Electronic Structure of Spin-Chain Compounds: Common Features 
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The incommensurate composite systems M14CU24O41 (M=Ca,Sr,La) are based on two fundamen- 
tal structural units: Cu02 chains and CU2O3 ladders. We present electronic structure calculations 
within density functional theory in order to address the interrelations between chains and ladders. 
The calculations account for the details of the crystal structure by means of a unit cell comprising 
10 chain and 7 ladder units. It turns out that chains and ladders can be treated independently, 
which allows us to introduce a model system based on a reduced unit cell. For the Cu02 chains, we 
find two characteristic bands at the Fermi energy. Tight binding fits yield nearest and next-nearest 
C"| , neighbour interactions of the same order of magnitude. 

a 

PACS numbers: 71.10.Pm: Fermions in reduced dimensions, 71.20.-b: Electron density of states and band 
£f*l , structure of crystalline solids, 74.72.-h: Cuprate superconductors 
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The isostructural spin-chain compounds M14CU24O41 (M=Ca,Sr,La) have been subject of intensive research in 
recent years, mainly due to their rich phase diagram and close relations to the high-T c cuprates. Their incommensurate 
crystal structures consist of planes of quasi one-dimensional CuC>2 chains stacked alternately with planes of two-leg 
CU2O3 ladders. The orientation of these chains and ladders defines the crystallographical c-axis, where the lattice 
constants of the two subsystems satisfy in a good approximation 10c c hain ~ 7ci a dder [3, 0j- The copper ions are 
intrinsically hole doped with nominal Cu valence +2.25 for both the Ca and the Sr compound. However, the ladders 
contain less holes than the chains. The resulting Cu valence of +2.75 for the chain sites corresponds to a quarter filled 
one-dimensional band [jj. Optical conductivity as well as x-ray absorption experiments suggest that substitution of 
Ca for Sr induces a transfer of holes from the chains back to the ladders [4| . On La substitution the intrinsic doping 
decreases, reaching the undoped state with nominal Cu valence +2.00 for La6CasCu2404i. 
O , In the case of Sri4Cu2404i the Cu02 chains are non-magnetic with a spin gap of about 130 K [5J and the CU2O3 
ladders show a charge density wave in a wide temperature range 0, 0|- In contrast, for Ca rich samples antifcrro- 
y—( ] magnetic ordering is found in the chains [H, 0| and the ladders become superconducting under pressure fiol ] . This 

■ strong doping dependence of both the charge and the magnetic order [ll|, [l2|, EH, 0] is accompanied by only minor 
modifications of the crystal structure (l5l flq. ll7T] . From the theoretical point of view, the magnetic properties of un- 
doped two-leg S=l/2 ladders have been analyzed be means of a Heisenberg model including cyclic ring exchange [l8l |. 
Possible effects of hole doping are investigated in A large spin gap makes the CU2O3 ladder planes magnetically 
inert. On the other hand, the magnetic phase diagram of the Cu02 chains is quite complicated. Ferromagnetism, as 

■ expected for Cu-O-Cu bond angles of approximately 90° in the chains, is realized in La rich systems. The origin of 
the observed intrachain antiferromagnetic order for a quarter filled Cu band, however, led to controversial discussions 

"j-^ ■ in the literature [13, H9> HH . Understanding of the doping dependence of both the charge order and the spin gap is 
likewise far from complete. 

Electronic structure calculations applying density functional theory and the local density approximation have been 
performed by Miiller et al. [22j for the isostructural CU2O3 ladders in SrCu203 and by Arai et al. [23j for a simplified 
Si'i4Cu2404i unit cell. Probably because of a high demand on CPU time, state-of-the-art band structure calculations 
accounting for the details of the crystal structures are absent so far. In the present letter we aim at filling this gap in 
\ order to address the interplay between the chain and ladder subsystems and to identify electronic features common 
to all compounds due to the specific crystal structure. Moreover, we introduce a model system capturing the relevant 
properties of the CUO2 chains for various compounds and levels of hole doping. By this model system, we identify 
the electronic states at the Fermi level and study intrachain nearest and next-nearest neighbour interactions in the 
framework of a tight binding model. As we account for the decisive structural ingredients, we obtain, as far as we 
know for the first time, realistic insights into the electronic structure of the Cu02 chains in spin-chain compounds. 

Our first principles band structure calculations for the full Sri4Cu2404i unit cell are based on the scalar-relativistic 
augmented spherical wave (ASW) method [24| , which has proven to be suitable for dealing with unit cells comprising 
a large number of atomic sites [2a, I2IIJ 27I]. As structural input we apply the crystallographical data due_to a recent 
structural refinement by Gotoh et al. 0- A previous investigation of the crystal structure by Ukei et al. [l[ seems to 
suffer from low quality samples. Since the Sri4Cu2404i unit cell contains two formula units, we have to deal with 28 
strontium, 48 copper, and 82 oxygen atomic spheres. Furthermore, 146 additional augmentation spheres are necessary 
in order to represent the correct shape of the crystal potential in the large voids of the open crystal structure. The 
basis sets taken into account in the secular matrix therefore comprise Sr 5s, 5p, Ad, Cu 3d, 4s, Ap, and O 2s, 2p 
orbitals, as well as states of the additional augmentation spheres. Brillouin zone integrations are performed using an 
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FIG. 1: Partial Cu 3d (first row) and O 2p (second row) densities of states (per atom) for the chain (left) and ladder (right) 
subsystems. The results are based on the full Sri4Cu2404i unit cell. 



increasing number of finally 72 k-points in the irreducible wedge. The Vosko-Wilk-Nusair parametrization is used for 
the exchange-correlation functional. 

The model system for the CuC>2 chain discussed subsequently comprises 6 atomic sites. It is studied via band 
structure calculations using the WIEN2k program package, which is a popular full-potential linearized augmented 
plane wave (LAPW) code [2g|. To obtain reliable results, we apply basis sets of 18135 plane waves and k-meshes with 
up to 108 points in the irreducible wedge of the Brillouin zone. Whereas the Cu 3p and O Is orbitals are treated 
as semi-core states, the valence states consist of Cu 3d, 4s, Ap and O 2s, 2p orbitals. The data are based on the 
Perdew-Burke-Ernzernhof scheme for the exchange-correlation functional. 

Figure Q] displays partial Cu 3d densities of states (DOS) resulting from our calculations for the full Sr^Ct^Cn 
unit cell. The Cu 3d DOS is seperated into contributions due to the Cu02 chains (left) and CU2O3 ladders (right), 
and normalized with respect to the number of contributing sites. We take into account the structural details of the 
compound by means of a unit cell comprising 10 chain and 7 ladder units along the crystallographical c-axis. This 
unit cell gives rise to an excellent approximation to the incommensurate crystal structure [2j. Note that different 
lengths of the chain and ladder units in principle lead to slight periodical tiltings, which do not affect our conclusions 
and therefore are not considered here. 

Both for the chains and ladders there is a finite DOS at the Fermi energy in fig. [H contradicting the experimentally 
observed non- metallic state of Sri4Cu2404i. In order to reproduce the insulating behaviour, it would be necessary 
to include electronic correlations beyond the local density approximation. However, for our purpose of studying the 
interplay between the chains and ladders such local interactions are not of interest. Electronic states in the energy 
range of fig.Q]are composed of Cu 3c? and O 2p orbitals (almost exclusively), which indicates strong Cu-0 hybridization 
both in the chains and in the ladders. As expected, this resembles the findings for the simplified unit cell of Arai et 
al. [HI]. The same is true for the gross shapes of the chain and ladder DOS curves. While the ladder DOS is spread 
over the whole energy range from — 7eV to 1.5 eV, the chain DOS shows a distinct band with a width of leV near the 



3 



O 
Q 



chains only 






o 


- 








6 - 












4 - 








w 






2 - 


) 










\ 












- 



-6 -5 -4 -3 -2 
(E - E F ) (eV) 



-1 1 



ladders only 





-7 -6 -5 -4 -3 -2-10 1 
(E - Bp) (eV) 



FIG. 2: Partial Cu 3d densities of states for chain (left) and ladder (right) copper sites. The results are based on Sri4Cu2404i 
unit cells where either the chain or the ladder substructure was removed. As concerns the vicinity of the Fermi level, the curves 
resemble almost perfectly the findings for the chain and ladder copper sites in the full Sri4Cu2404i unit cell, respectively, 
compare fig. [T] 
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FIG. 3: Partial Cu 3d densities of states for a model Cu02 chain (left) and two coupled Cu02 chains (right), arranged 
symmetrically and unsymmetrically. The states are subject to a rigid band shift of about 1.5 eV as compared to the real 
system. 



Fermi level, separated by a gap of likewise 1 eV from the other valence bands. This band is close to quarter filling, as 
to be expected, and hence subject to a large variety of possible ordering processes of the charge and spin degrees of 
freedom [2^ |. 

Turning to the question about the interrelations between the chain and ladder subsystem, the Sr^C^,^^ unit 
cell is modified as follows: (a) We remove either the chain or the ladder substructure and (b) replace Sr by K. Since 
the M ions in M14CU24O41 act almost exclusively as charge donors, the latter allows us to maintain the correct doping 
when removing one of the substructures. For the full unit cell, formally one electron is transferred from each Sr site to 
both the chain and the ladder subsystem. The same amount of charge transfer is achieved for the half-filled unit cell 
via the potassium, because hybridization effects between Sr/K and Cu/O orbitals are negligible. Figure [5] displays 
DOS curves for the chains only configuration on the left and for the ladders only configuration on the right side. In 
both cases the findings resemble the corresponding curves of fig. [21 as concerns the total width of the valence bands 
as well as the gross shape of the DOS. In the vicinity of the Fermi energy even the details of the DOS are identical. 
In particular, for the Cu02 chains, the band width and the filling of the distinct band around the Fermi level hardly 
alter. 

To conclude, since the above argumentation does not depend on the specific choice of the donor ion M, we find 
that spin-chain compounds M14CU24O41 are well described in terms of three largely independent subsystems: The 
Cu02 chains, the CU2O3 ladders, and the electron donor system of M ions, which separates chains and ladders. 
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FIG. 4: Electronic band structure (Cu 3d bands) for coupled CUO2 chains, arranged symmetrically (left) and unsymmetrically 
(right). The results of tight-binding fits are summarized in tabel|T] 





band 1, symm. 


band 2, symm. 


band 1, unsymm. 


band 2, unsymm. 


ti (eV) 


-0.115 


-0.105 


-0.135 


-0.080 


t 2 (eV) 


-0.183 


-0.098 


-0.165 


-0.115 


eo (eV) 


1.485 


1.645 


1.460 


1.675 



TABLE I: Tight-binding parameters for the Cu 3d bands of fig. [4] 



Hybridization between atomic orbitals belonging to different subsystems can be neglected because the valence states 
of the M ions are almost fully depopulated. As a consequence, we attribute charge transfer between chains and ladders 
when replacing Sr by isoelectronic Ca [H, Q to simultaneous modifications of the crystal structure. 

Since we have established that the structural subsystems of spin-chain compounds can be treated independently, we 
address the electronic features of the CuC>2 chains in the following via a model system based on a reduced unit cell. 
For convenience, we choose the Euklidean z-axis as the chain axis and use the interatomic distances and bond angles 
of Sri4Cu2404i. Each Cu site thus is connected to two symmetrical O sites in both chain directions, where the bond 
lengths are 1.89A in one and 1.87 A in the other direction. The related O-Cu-0 bond angles are 85.5° and 86.7°, 
respectively. Neighbouring Cu atoms are separated by 2.75 A. In order to prepare for the symmetry analysis of the 
Cu bands, we align the CuC>2 chains in the Euklidean a;z-plane. As a consequence, Cu 3d xz orbitals mediate the main 
part of the Cu-0 overlap. We study both a single Cu02 chain, for which the unit cell comprises just one CuC>2 unit, 
and two coupled chains. The latter allows us to address the influence of relative shifts between adjacent chains with 
respect to the chain axis, which distinguish specific spin-chain compounds. For Cai3.gSro.4Cu2404i [9I [Til] the chains 
are shifted by half their intrachain Cu-Cu distance, for example, while the shift is only 30% of the Cu-Cu distance for 
Sri4Cu2404i. We compare these two cases, which we call the symmetrically and unsymmetrically arranged chains, 
respectively. 

Figure [3] summarizes the results of our LAPW model calculations, showing partial Cu 3d densities of states for a 
single Cu02 chain on the left side and two coupled CuC>2 chains on the right side. All three DOS curves resemble the 
essential features of the ASW data for the real system, see figs. Q] and [2] We obtain a widespread structure with a 
width of about 5.2 eV at lower energies and a structure extending over a range of leV at higher energies. However, 
comparison with the ASW findings results in a rigid band shift of slightly less than 1.5 eV to higher energies. This shift 
traces back to the fact that the M ions are not taken into account in the model system, see the previous discussion. 
In the following we hence are interested in the electronic states in the energy range from about 1 eV to 2 eV. As 
expected, these states are mainly due to the antibonding combination of Cu 3d xz and both O 2p x and 2p z orbitals. 
Contributions of other states amount to less than 0.1% of the total DOS at the Fermi energy. In order to obtain 
further insight into the electronic features of the Cu02 chains, we turn to the band structure underlying the DOS 
curves of fig. [3] 

Because the Cu02 chains run parallel to the z-axis of our unit cell, fig. 2] depicts the Cu d xz bands in the corre- 
sponding direction of the orthorhombic Brillouin zone, i.e. between the high symmetry points T and Z. As the length 
of a Cu02 unit amounts to c c hai n =2.75 A, we can write kr = (0,0,0) and k^ = (0, 0, 27r/c c h a in)- Reasonable fitting of 
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the bands is possible by assuming a the tight-binding dispersion of the form 

e(fcz) = eo - 2£i cos(fc 2 c chain ) - 2i 2 cos(2fc z c chain ), (1) 

where k z is the ^-component of the reciprocal lattice vector and eo is the band center. Moreover, t\ and t% denote 
the nearest and next-nearest neighbour hopping parameters along the chains, respectively. The outcome of the 
tight-binding fits for the bands of fig. [4] is given in table Q] Beyond our one-dimensional model (JTJ) , influence of finite 
interchain hopping can be identified in the band structure for the symmetrical configuration, but is not relevant for our 
further conclusions. In contrast to Arai et al. [23] , we find two Cu 3d xz bands, corresponding to the two coupled chains. 
We attribute this difference to inadequate crystallographical data entering the former band calculation. Importantly, 
we obtain for both Cu 3d xz bands nearest and next-nearest neighbour couplings of the same order of magnitude, 
independent of shifts between neighbouring CuC>2 chains. Despite the 90° Cu-O-Cu intrachain bond angles, strong 
nearest neighbour interaction therefore appears to be typical of the whole class of spin-chain compounds. As the 
calculated magnitudes of the nearest and next-nearest neighbour coupling are common to all systems, we conclude 
that changes in the magnetic ordering on replacing Ca by La result almost exclusively from the modified filling of the 
Cu 3d xz bands. 

In summary, we have studied the electronic structure of Sri4Cu2404i by means of density functional theory. Taking 
into account the details of the crystal structure, we have addressed the interrelations between the structural subsystems 
the spin-chain compounds M14CU24O41 (M=Ca,Sr,La) are composed of. Hybridization between the CuC>2 chains, the 
Cu 2 03 ladders, and the electron donor M ions is found to be negligible. Since these subsystems consequently can be 
treated individually, we have analyzed a simplified model system for the Cu02 chains to identify features common 
to all spin-chain compounds. It turns out that two characteristic Cu 3d xz bands at the Fermi energy dominate the 
electronic and magnetic properties of the chains. The filling of these bands distinguishes different compounds. Tight- 
binding fits prove that nearest and next-nearest neighbour interactions in the Cu02 chains are of the same order of 
magnitude, which applies to the whole class of spin-chain systems. 
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